Relaxation in statistical many-agent economy models 
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We review some statistical many-agent models of economic and social systems inspired by mi- 
croscopic molecular models and discuss their stochastic interpretation. We apply these models to 
wealth exchange in economics and study how the relaxation process depends on the parameters of 
the system, in particular on the saving propensities that define and diversify the agent profiles. 
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I. INTRODUCTION 



One might question how theories that try to explain 
the physical world of elementary particles, atoms, and 
molecules can be applied to understand the social struc- 
ture in its complexity and the economic behavior of hu- 
man beings: is it possible to describe the behavior of 
people with simple models? Is it even possible to iden- 
tify and quantify the nature of the interactions between 
them? Even though it is still difficult to find answers to 
these questions, during the past decade physicists have 
made attempts to study problems related to economics, 
the social science that seeks to analyze and describe the 
production, distribution, and consumption of wealth [1]. 

Here we will not try to review all these attempts, rather 
we briefly describe what we name statistical many-agent 
models. In these models, economic activity is described 
as a flow of wealth between basic units, referred to as 
agents, representing e.g. individuals or companies. Each 
of the N agents {1, 2, . . . , i, . . . , N} has a wealth 
Xi, that changes in time as agents exchange wealth be- 
tween each other, according to the trading rules detailed 
in Sec. II. These underlying trading rules only depend 
on one set of parameters, namely the saving propensities 
{Xi}, with < A< < 1. Statistical many- agent models de- 
scribe closed economy systems and can reproduce some 
features of wealth distributions, such as an exponential 
at intermediate values of wealth and a power law at high 
values. In a particular model, all agents have the same 
global A; in a more general model, different agents have 
different values of A. 

If the saving propensity is equal for all the units 
(global saving propensity models), Xi = A, the equilib- 
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rium wealth density f(x) is given by a T-distribution, 



f(x) = fl ln {Px) 



JL 
r(n) 



D 1 + 2A 



1 - A 



{pxT' 1 exp(-/3x) , (1) 
(2) 



Here n is a real number in the interval [1, oo) and D = 2n 
can be considered as the effective dimension of the sys- 
tem: in fact the distribution j n ([3x) is just the Maxwell- 
Boltzmann distribution for the kinetic energy x of a gas 
in thermal equilibrium at a temperature /3 _1 in a D- 
dimensional space. Thus the parameter /? can be inter- 
preted as the inverse temperature: consistently with the 
equipartition theorem, /3 _1 — 2(x)/D = (x)/n, where 
the constant (x) is the average wealth (x) = X/N and 

X = 53j=i x i is the total wealth. In an economic system, 
temperature is proportional to the fluctuations of wealth 
around its average value. The model with a single global 
saving propensity describes well wealth distributions at 
intermediate values of x, but cannot predict power laws 
at large x. 

If A is uniformly distributed among the agents accord- 
ing to a given density </>(A) (distributed saving propensity 
models), then one finds, under quite general conditions 
on the shape of <^(A), an exponential law at intermediate 
values of x and a robust Pareto law, 



f(x) oc x' 



(3) 



with a > 1, at large x. This power law was suggested by 
Vilfrcdo Pareto [2] more than a century ago to describe 
the tail of wealth distributions and is usually found to be 
characterized across various countries by a Pareto expo- 
nent a ps 3/2. 

The basic exchange laws underlying some of these mod- 
els are reviewed in Sec. II. In Sec. Ill we focus on the 
relaxation to equilibrium generated by the exchange laws 
of models with uniformly distributed A and illustrate its 
dependence on the saving propensity A through some ex- 
amples. We also consider the corresponding relaxation 
time distributions and discuss the relation between the 
Pareto exponent and the relaxation behavior of the sys- 
tem. Results are summarized in Sec. IV. 
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II. MANY-AGENT MODELS OF A CLOSED 
PURE EXCHANGE ECONOMY 

Our aim is to study a general many-agent statistical 
model of a closed economy without growth (analogous to 
the kinetic theory model of ideal gases, see Fig. 1), where 
N agents exchange a quantity x, that we have defined as 
wealth. The states of the agents are specified only by 
the wealths {x{\, while the total wealth X = Yli=i x i 
is conserved. The evolution of the system is carried out 
according to the following algorithm: at every time step 
two agents i and j are chosen randomly and an amount 
of wealth Aa; is exchanged, so that the agent wealths x\ 
and x'j after the transaction are 

■2/^ — Xi . 

x'j = Xj + Ax . (4) 

Different transaction rules have been studied analytically 
or numerically by various authors and are summarized 
here below. 



A. Basic model without saving: Boltzmann 
distribution 

A stochastic trading rule, that redistributes the 
wealths of two agents randomly, was introduced in Ref. 3, 

Xj ^(^2 ~h : 

x'j = l(xi +Xj) , (5) 

where e is a uniform random number in (0, 1) and e + e = 
1. Eqs. (5) are equivalent to the trading rule of Eq. (4) 
if Ax — exi — exj. In another version of the model, 
the money difference Ax is assumed to have a constant 
value independent of the two trading agents [4-6], i.e. 
Aa; = Axo. Both these forms for Aa; lead to a robust 
equilibrium Boltzmann (or Gibbs) distribution, 

f{x) = /T 1 exp(-#r) , (6) 

with effective temperature = (x) = X/N [3-6]. 

B. Model with a global saving propensity: Gamma 
distribution 

Introducing a saving criterion through a saving propen- 
sity parameter < A < 1 [7, 8] modifies the trading rule 
as follows: 

x\ — Xxi + e(l - \)(xi + xj) , 

x'j — Xxj + e(l — X)(xi + Xj) . (7) 

This model is similar to that introduced by John Angle 
in 1983 on the basis of the Surplus Theory [9-11], which 
however differs both in the mathematical definition of the 




FIG. 1: Analogy between the minimal closed pure exchange 
economic model and a classical isolated system of an ideal gas. 
In the latter case particles undergo random elastic collisions 
and exchange a fraction of their kinetic energy, while in the 
closed economy model agents perform random trades with 
each other and exchange a fraction of their wealth according 
to some statistical rule. 



exchange rule and interpretation [11]. A closer compari- 
son will be studied elsewhere, while here, for the sake of 
simplicity, we will focus on the trading rule (7) and some 
of its generalizations for studying the relaxation process. 
The rule of Eq. (7) corresponds to the process defined by 
Eq. (4) if 

Ax = (1 - X)(eXi - exj) . (8) 

The parameter A represents the fraction of wealth saved 
before the reshuffling takes place. The resulting equilib- 
rium distribution is qualitatively different from a simple 
exponential function, being a T-distribution [12, 13], see 
Eq. (1), that has a mode x m > and a zero limit for 
x — > 0, see Fig. 2. 

C. Models with a continuous distribution of saving 
propensities 

More realistic and interesting models are obtained 
when agents i = 1, . . . , N are diversified by assigning 
them different saving propensities Aj [14-21], e.g. with 
the Xi distributed uniformly on the interval [0,1). The 
trading rule is then 

x\ = XiXi + e[(l - Xi)xi + (1 - Xj)xj] , 

x'j = XjXj + e[(l - A i )a; l + (1 - Xj)xj] , (9) 

or, equivalently, can be formulated through Eq. (4) with 
A a; given by 

Ax = e(l - Xi)xi - e(l - Xj)xj . (10) 

Numerical simulations and theoretical considerations 
suggest that these models relax toward a robust power 
law oc l/a; 1+Q with a Pareto exponent a = 1 in the case 
of uniformly distributed A and with a > 1 if the A-density 
0(A) ->0as 0(A) - (1 - A)"" 1 for A — > 1. In the fol- 
lowing we study the relaxation process of models with 
uniformly distributed A. 
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FIG. 2: Wealth probability density in linear (left) and semi-log (right) scale for various global saving propensities A. The 
exponential curve A = is the equilibrium solution for the basic model (Sec. II A), while the other curves correspond to a 
nonzero global saving propensity (Sec. II B). 



III. RELAXATION PROCESS 

If a real economic system is characterized by a wealth 
distribution with a certain shape, it is of great interest to 
know on which time scale the system relaxes toward this 
distribution from a given arbitrary initial distribution of 
wealth, and how the relaxation process depends on the 
system parameters, in particular on the system size and 
the distribution of saving propensities. 

In the simulations presented below, all agents start 
from the same initial wealth Xi(t = 0) = xq = 1. The 
value xq, due to the conservation of the total wealth 
X = JDtLi Xi, also represents the global average value 
of x at any time t, i.e. (x(t)} = J xf(x) dx = x = X/N. 
This setup is used to model a more general situation 
where the initial conditions of the agents are far from 
equilibrium. 
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FIG. 3: Average wealth {x{t))\-o.Q9 versus the rescaled time 
t — T/N for systems with different number of agents N = 
100, 1000, and 10 000, but the same saving propensity density 
cA(A), that uniformly partitions agents into 100 subsets with 
A = 0.00, 0.01, . . . , 0.99. T is the total number of trades. 



A. Relaxation to equilibrium as a function of 
system size 

Before analyzing the dependence of the time scale on 
the saving propensity distribution, we shortly consider 
its dependence on the number of agents N. If time is 
measured by the number of transactions T, we find that 
the time scale is proportional to the number of agents 
N: a system A that is m times larger than a system 
B (A^a = itiNb) relaxes m times slower than B. This 
is shown in Fig. 3, where the average wealth (x(t))\ of 
the agent subset with A = 0.99 is plotted for various 
systems with different values of N versus the rescaled 
time t = T/N. However, the A-density </>(A) is the same 
for all systems and uniformly partitions each system into 
100 subsets with values A = 0.00, 0.01, . . . , 0.99. 

Here and in the following we define time t as the ratio 
t = T/N between the total number of trades T and the 
total number of agents N, i.e. what is usually called a 
Monte Carlo cycle or sweep in molecular simulation [22]: 



in a Monte Carlo cycle, each agent performs on average 
the same number of trades (actually two), in the same 
fashion as in molecular dynamics each particle is moved 
once at every time step. The results will not change if 
one of the two agents involved in an exchange is selected 
sequentially, e.g. in the order of its index i = 1, . . . , N, 
as is common practice in molecular simulations. This en- 
sures that every agent performs at least one trade per 
cycle and reduces the amount of random numbers to be 
drawn. The previous considerations suggest the intro- 
duction of a time unit To, such that during any time in- 
terval (i, t + tq) all agents perform on average one trade 
(or the same number of trades). In this way the dynam- 
ics and the relaxation process become independent of N. 
The existence of a natural time scale independent of the 
system size provides a foundation for using simulations 
of systems with finite N in order to infer properties of 
systems with continuous saving propensity distributions 



4 



and N — > oo. 



B. Relaxation to equilibrium as a function of 
saving propensity 

Relaxation in systems with constant A has already 
been studied in Ref. 7, where a systematic increase of the 
relaxation time with A, and eventually a divergence for 
A — > 1, was found: for A = 1, no exchanges can occur, so 
that the system is frozen. Here we consider systems with 
uniformly distributed A. In this case a similar behavior 
of the relaxation times is observed, broken down to sub- 
systems with similar values of A. As discussed in detail in 
Refs. 21, 23, 24, the partial wealth distributions of agents 
with a given value of A relax toward different states with 
characteristic shapes f\(x). The generic function f\(x) 
has a maximum and an exponential tail, thus closely re- 
calling the shape of a T-distribution. The corresponding 
average value is given by (x)\ = J xf\(x) dx = fc/(l — A), 
where k is a suitable constant determined through the 
condition J (x) \ 0(A) dX = X/N; X is the total wealth 
of the system. Even if the partial distributions decay 
exponentially with x, the sum of all partial distribu- 
tions results in a Pareto law at large values of x, i.e. 
f( x ) = J2\f*( x ) ~ l/x 1+a . Numerical simulations 
clearly show that agents with different values of A are 
associated to different relaxation times t\. 
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FIG. 4: Mean wealth {x(t))\ versus time for various A. The 
10 4 agents are uniformly partitioned into 100 subsets with 
A = 0.01, 0.02, . . . , 0.99. Higher A correspond to longer 
relaxation times. The continuous line (x) — xq partitions 
agents into poor (x < xq) and rich (x > xq) ones. 



Results are illustrated in Fig. 4 for a system of N = 
10 4 agents uniformly partitioned into 100 subsets with 
A = 0.01, 0.02, . . . , 0.99: mean wealths of subsets cor- 
responding to a value of A closer to 1 relax slower toward 
their asymptotic average wealth (x}\ oc 1/(1 — A). 

The average wealth xq allows to introduce a thresh- 
old that partitions the system into poor agents, with an 
asymptotic average wealth {x(t — > oo))a < xo, and rich 
agents with (x(t — > oo))a > Xq. The poor-rich threshold 



(x) = xo — 1 is represented as a continuous line in Fig. 4 
and corresponds to A rj 0.75 for this particular example. 

The differences in the relaxation process can be re- 
lated to the different relative wealth exchange rates, that 
by direct inspection of Eqs. (9) and (10) appear to be 
proportional to 1 — A. Thus, in general, higher saving 
propensities are expected to be associated to slower re- 
laxation processes. A more detailed analysis can be car- 
ried out as shown in Fig. 5: after the rescaling of time and 
wealth by the factor (1 — A), mean wealths corresponding 
to agents with different values of A (Fig. 5 left) appear to 
relax approximately on the same time scale and toward 
the same asymptotic value (Fig. 5 right). In fact, the fac- 
tor (1 — A) is proportional to the wealth exchange rates 
and, at the same time, through the condition of station- 
arity, determines the equilibrium average wealth values 
(x)\ = k/(l — A) [21]. Agents start from the same ini- 
tial condition Xi{t = 0) = xq = 1. In this case, in order 
to study in greater detail the high saving propensity pa- 
rameter region, that corresponds to the high relaxation 
time region, the system of N = 10 4 agents has been uni- 
formly partitioned into 200 subsets with saving propen- 
sities A = 0.5000, 0.5025, . . . , 0.9975. Actually, this is 
not a uniform distribution of A on [0, 1), since 0(A) = 
for A < 0.5, however it does not matter because what 
counts is the high saving propensity parameter interval. 

C. Relaxation time distribution 

The model with distributed saving propensities is com- 
pletely specified by the trading rules of Eqs. (9) and the 
set of saving propensities {A^} of the N agents. In the 
case of a continuously distributed A, a continuous sav- 
ing propensity density 0(A) can be used in place of the 
discrete A-set, normalized so that J Q 0(A) dA = 1. 

Here we suggest a method to obtain the wealth as well 
as the relaxation distribution directly from the saving 
propensity density 0(A). It follows from probability con- 
servation that f(x)dx = 0(A)dA, where a; is a short no- 
tation for (x)\ and f{x) is the density of the average 
wealth values. In the case of uniformly distributed sav- 
ing propensities, one obtains 

/» = *<f -*('-!)?■ <") 

This shows that a uniform saving propensity distribution 
leads to a power law f(x) ~ 1/x 2 in the (average) wealth 
distribution. In general a A-density going to zero for 
A -> 1 as 0(A) oc (1 - A) a_1 (with a > 1) leads to the 
Pareto law f(x) ~ l/x 1+a with Pareto exponent a > 1 
as found in real distributions. 

In a very similar way it is possible to obtain straight- 
forwardly the associated distribution of relaxation times 
■0(t) for the global relaxation process through the rela- 
tion between the relaxation time t\ and the agent saving 
propensity: given that the time scale follows a relation 
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FIG. 5: Left: average wealth (x(t)}\ versus (log of) time for the A's listed in the figure. The 10 4 agents are partitioned into 200 
subsets with A = 0.5000, 0.5025, . . . , 0.9975. Right: same as in the left figure after rescaling wealth and time by the factor 
(1 — A), inversely proportional to the mean wealth and proportional to the average wealth exchange rate. 



t\ oc 1/(1 — A), then 

where r' is a proportionality factor. Comparison with 
Eq. (11) shows that and f(x) are characterized by 
power law tails in r and x respectively with the same 
Pareto exponent. 

It is to be noticed, as discussed in Ref. 21, that in the 
parameter region A — > 1, from which the main contri- 
butions to the Pareto power law tail come, the widths 
of the generic equilibrium partial distributions f\(x) in- 
crease more slowly than the difference between the mean 
values (x)\> — (x)\ corresponding to two agents with con- 
secutive values of the saving propensity A' and A. This 
implies that at equilibrium and in the tail of the distribu- 
tion it is possible to resolve the mixture J2\ fx( x ) m t° its 
components f\(x) and to approximate the current value 
of wealth x(t) of a certain agent with saving propensity 
A (that is actually a stochastic process) with the corre- 
sponding average value, (x)\ « x, so that f(x) ~ f{x). 

Finally we notice that an ensemble with a power law 
distribution of relaxation times undergoes a slow relax- 
ation process if the exponent of the relaxation time dis- 
tribution is smaller than two, so that a Pareto exponent 



larger than two, as automatically generated by the model, 
seems to ensure a normal relaxation. 



IV. CONCLUSIONS 

The relaxation process of statistical many-agent mod- 
els of a closed pure exchange economy, where trading is 
described as a flux of wealth between different agents, has 
been found to be slower for agents with a larger saving 
propensity parameter A, who are also the agents resulting 
richer at equilibrium. For a uniform A-distribution, the 
relaxation time is t\ oc 1/(1 — A). Furthermore, a smooth 
distribution of saving propensities leads to distributions 
of wealth and of relaxation times characterized by power 
law tails with the same Pareto exponent a > 1, which 
ensures a fast relaxation toward equilibrium. 

We also remark that if time is measured in Monte Carlo 
cycles, i.e. the ratio between the total number of trades 
and the total number of agents, so that every agent per- 
forms on average two trades during a cycle, the time 
evolution and the relaxation process are independent of 
the system size, thus providing information on arbitrarily 
large systems. 
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